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In view of promising applications of fractal nanostructures, we analyze the spectra of quantum 
particles in the Sierpinski carpet and study the non-correlated electron gas in this geometry. We 
show that the spectrum exhibits scale invariance with almost arbitrary spacing between energy 
levels, including large energy gaps at high energies. These features disappear in the analogous 
random fractal—where Anderson localization dominates—and in the regular lattice of equally sized 
holes—where only two length scales are present. The fractal structure amplifies microscopic effects, 
resulting in the presence of quantum behavior of the electron gas even at high temperatures. Our 
results demonstrate the potential of fractal nanostructures to improve the light-matter interaction 
at any frequency, with possible applications, e.g., in the development of solar cells with a wide 
absorption spectrum, artificial photosynthesis, or nanometamaterials with tailored Fermi levels and 
band gaps, operating in a wide range of frequencies, and with extended operating temperature range. 


New generation of materials, based on nanostructures, 
is showing its potential to revolutionise industry [1, 2]. 
In parallel, self-similarity properties of fractal structures 
have revised the way we manufacture antennas [3] and 
permitted the production of novel metamaterials [4, 5]. 
New materials and manufacturing techniques have en¬ 
abled producing fractal circuits [6], porous media such as 
silicas or aerogels, where superfluid helium flows in a frac¬ 
tal environment [7], fractal-shaped clusters of nanoparti¬ 
cles formed by aggregation [8]; or even fractal snowflakes 
composed of graphene [9]. Additionally, one can find self¬ 
similarity in the electronic spectra of quantum systems, 
such as in Hofstadter’s butterfly and in the fractal quan¬ 
tum Hall effect [10, 11]. Indeed, quantum fractals are 
drawing an increasing attention due to their intriguing 
properties and potential applications [12]. 

To describe the properties of a quantum particle in a 
fractal geometry (e.g., an electron trapped in a fractal¬ 
shaped nanowire), we first find the eigenstates of the par¬ 
ticle’s Hamiltonian 77 = T + V, where T is the kinetic 
energy operator and V is the potential energy of an in¬ 
finite well with the desired fractal boundaries. For the 
well, we have chosen the paradigmatic Sierpinski carpet, 
a fractal object studied in many different contexts [6, 13]. 
This carpet is constructed by subdividing a square into 
smaller copies of itself and by removing the central copy 
iteratively, as shown in Fig. 1. Each iteration is labeled 
with an integer G indicating the generation of the carpet. 
We study generations up to G = 6, starting from G = 0, 
which represents the square well. Quantum Hamiltonian 
for a particle of mass m in a fractal geometry can be 
written as PL = T + V = + V(q), where V(q) is 

an infinite well potential with the shape of the accessible 
region in the carpet of generation G and side L. With¬ 
out loss of generality, units of L for length, K 2 /mL 2 for 
energy and ksh 2 /mL 2 for temperature (where ks is the 
Boltzmann constant) are used when m, H, or L are not 
shown explicitly. For any finite generation G, Sierpinski 


carpet can be viewed as a polygonal billiard with rational 
angles, and therefore is an example of pseudointegrable 
systems discovered by Richens and Berry [14], who also 
studied the equivalent of the first generation of Sierpin¬ 
ski carpet, albeit with periodic boundary conditions. As 
for higher fractal generations, Katomeris and Evangelou 
studied properties of Sierpinski lattices [12]; we instead 
focus on the dynamics of quantum particles in a frac¬ 
tal Sierpinski volume, such as electrons in fractal shaped 
wires, which are not yet well understood. 

To obtain the eigenenergies and eigenstates of Li, we 
have used two different diagonalization methods: the 
finite differences method (FDM) and imaginary-time 
nonuniform-mesh method (ITNUMM) [15]—see section I 
in Supplementary Material (SM) for details. The relative 
error in the eigenenergies computed with FDM grows lin¬ 
early with the energy, whereas the error remains roughly 
constant for the ITNUMM. Thus, despite the error of 
the FDM being smaller than the error of ITNUMM for 
the lowest-energy states, it eventually becomes larger for 
high-energy states. As a consequence, ITNUMM is the 
preferred method if one is interested in energies of many 
states. Nevertheless, we have verified all our numerical 
results by evaluating the energies both with ITNUMM 
and FDM. 

The spinless single-particle spectrum. Figure 1 shows 
the structure of the lowest-lying states from the zeroth 
to fourth generation (G = 0, 1, 2, 3 and 4. Gener¬ 
ations G = 5 and 6 are not shown due to the limi¬ 
tations of image resolution.) As the fractal generation 
increases, the wavefunctions develop intricate structures 
on smaller and smaller scales, resulting in the exponen¬ 
tial increase of the kinetic energy of the ground state 
e gs (Fig. 2a). This contribution to the energy affects 
equally also all excited states and manifests itself as a 
global shift of the energy spectrum for a given genera¬ 
tion. Hence, the mode number function A f (Fig. 2b) and 
energy levels e (Fig. 2c) are comparable for all genera- 
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Figure 1. Top panel: Iterative procedure for constructing each 
generation G of the Sierpinski carpet (“Fractal”), and of the 
analogous random fractal (“Random”) and regular non-fractal 
(“Regular”) structures. Bottom panel: Wavefunctions of the 
first five states (rows n = 1, 2, 4, 5) for the particle in a box 
(G = 0), and for the first four generations of the Sierpinski 
carpet (G = 1, 2, 3, 4). For comparison, corresponding results 
are shown for the third generation (G = 3) of both the random 
fractal carpet (where Anderson localization takes place) and 
the regular non-fractal lattice. 

tions when referred to the ground state. For G = 0, one 
obtains the spectrum of a particle in a two-dimensional 
well [16], where the mode number is proportional to the 
energy as A/"(e) = e mL 2 /2nli 2 , as predicted by Weyl’s 
law. Increasing the fractal generation gives rise to two 
qualitatively different energy regimes: (i) At lower ener¬ 
gies, some levels attract each other, agglomerating into 
bands, while others exhibit repulsion, generating large 
gaps between the bands; (ii) at higher energies, the 
mode number eventually approaches the prediction of 
Weyl’s law Af(e) = e m* L 2 /2 ttH 2 , where m* = (8/9 ) G m 
is the effective mass. Using this analytical form to 
rescale mode numbers as J\f* = J\f / 8° and energies as 
e* = emL 2 /9 G 2nh 2 (see Fig. 2d, left panel), we find, re¬ 
markably, that for any generation G the bands and gaps 
in regime (i) appear in the part of the spectrum with 
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Figure 2. Spectrum of a single particle in a Sierpinski-shaped 
well, a, Ground state energy e gs : The exponential decrease of 
the smallest length scale with increasing generation G mani¬ 
fests itself in the exponential increase of e ga . b, Mode number 
function Af(e — e ga ). c, Energy level diagram: The border be¬ 
tween low-energy, band-gap regime (black line segments) and 
the almost continuous high-energy regime (orange line seg¬ 
ments) satisfying Weyl’s law (see text) is shown for each G 
by a thick red line segment, d, The self-similarity of the 
spectrum becomes apparent when the energies and levels are 
rescaled according to Weyl’s law (left panel: Sierpinski carpet; 
right panel: random carpet). 


energies e* < 1 (AT* < 1), whereas the approximately 
constant density of states in regime (ii) appears for en¬ 
ergies e* > 1 (A f* >1)- This scaling of the energy levels 
demonstrates that the self-similarity in the geometry is 
reflected also in a scale-invariant spectrum. 

To further explore this phenomenon, we studied the 
eigenvalue nearest-neighbor spacing distribution in both 
regimes for all six generations. This distribution is 
widely used in the fields of random matrix theory and 
quantum chaos for characterizing ensembles of random 
matrices and dynamical properties of Hamiltonian sys- 
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terns [17, 18]. In addition, since the energy of any optical 
transition of an electron can be described as a sum of 
level spacings, their distribution can provide some infor¬ 
mation about the expected absorption/emission proper¬ 
ties of fractal nanowires. After first unfolding the energy 
levels (as described in section II of the SM), we computed 
the spacings s of every two consecutive modes as well as 
the maximum likelihood estimates of the parameters a 
and 6 of the Brody distribution [19] 

p(s) = (1 + b)as b exp (—as b+1 ) (1) 

fit to the probability density of s. Whereas Brody distri¬ 
bution with zero Brody parameter (6 = 0) reduces to the 
exponential distribution p(s ) = a exp (—as), reflecting 
Poisson level statistics characteristic of generic integrable 
systems [20], unit Brody parameter (b = 1) yields the 
Wigner distribution p(s) = 2asexp (—as 2 ) and, accord¬ 
ing to the quantum-chaos conjecture [21-24], reflects clas¬ 
sically chaotic dynamics of systems with time-reversal in¬ 
variant Hamiltonian. Generic time-reversal invariant sys¬ 
tems with mixed dynamics and without symmetry can be 
well described with a mixture of exponential and Wigner 
distributions [24, 25]; while this mixture itself is not a 
Brody distribution, it can be typically well fitted with 
a Brody distribution with 0 < b < 1. Remarkably, as 
shown in Fig. 3, for the Sierpinski carpet we find that 
the Brody parameter evolves from b = —0.03 ± 0.07 for 
G = 1, corresponding to Poisson level statistics (6 = 0), 
to 6 = —0.96±0.03 for G > 6 in regime (i), resulting in al¬ 
most an inverse power law distribution (6 = — 1). As the 
generation increases, the level repulsion-attraction inten¬ 
sifies and it becomes harder to define the mean level spac¬ 
ing. Indeed, for the asymptotic value 6 = — 1 the mean 
level spacing ceases to exist and the distribution cannot 
even be normalized; spacings of any order of magnitude 
are equally probable at any energy. This transition does 
not appear in regime (ii), where Brody parameter 0 
in all the cases studied here. While transitions between 
non-negative Brody parameters (from 6 = 1 to 6 = 0) 
have been seen in the family of Koch fractals of chang¬ 
ing fractal dimension [26], to the best of our knowledge a 
transition from Poisson statistics to a long-tailed distri¬ 
bution, similar to the one described here, was previously 
observed only by Katomeris and Evangelou, who saw a 
power law distribution with exponent —1.56 in a frac¬ 
tal lattice by removing diagonal disorder [12]. However, 
negative Brody parameters were observed previously in 
superintegrable systems (6 « —0.24) [27] or ID lattices 
with unbounded quantum diffusion (6 = —3/2) [28]. For 
the desymmetrized system (see SM section II), we find a 
qualitatively similar behavior involving a transition from 
positive to negative values, confirming that the effect is 
due to the fractal geometry and not due to degeneracies 
from simple geometrical symmetries. 

This peculiar feature of the spectrum of the Sierpin¬ 
ski carpet disappears (Fig. 3 and SM Fig. 9) both in the 
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Figure 3. a, dependence of the Brody parameter b of the dis¬ 
tribution of level spacings on the fractal generation G. Solid 
black: the low-energy regime (i) (see text) undergoes a transi¬ 
tion from the exponential (b, for G = 1, b ~ 0 reflecting Pois¬ 
son level statistics) to inverse power-law distribution (c, for 
G > 6, b « — 1); Solid gray: the high-energy regime (ii) shows 
no (or at most weak) dependence on G, approximately fol¬ 
lowing the Poisson distribution for all G. Dashed black: low 
energy regime (i) for the random fractal, evolving from the 
Wigner (6 ss 1) to exponential (b ta 0) distribution; Dashed 
gray: regime (ii) showing no (or at most weak) dependence 
on G > 2. 


analogous random fractal carpet, where the holes are dis¬ 
tributed randomly without overlapping (Fig. 1), and in 
the equivalent regular but non-fractal geometry, where 
all the holes are of the same size and are distributed reg¬ 
ularly in a lattice (Fig. 1). 

The random fractal exhibits a quasi-continuous spec¬ 
trum instead of the level repulsion and attraction pattern. 
Even if the potential well exhibits scale invariance in a 
statistical fashion, breaking the geometrical symmetries 
results in Anderson localization [29], typical of disordered 
media and here reflected in the localized eigenfunctions 
in Fig. 1. After rescaling the mode numbers and ener¬ 
gies as for the fractal, we find the same regime at high 
energies e* > 1 ( J\f * > 1, Fig. 2d right panel). However, 
in this regime for G = 2 and 3 the Brody parameter is 
6 = 0.90 ± 0.06 and b = 0.83 ± 0.12 respectively (Fig. 3), 
close to the value 1 characteristic of chaotic systems. On 
the other hand, the low-energy part of the spectrum— 
where Anderson states are found—lacks the band and 
gap structure, and the Brody parameter evolves, remark¬ 
ably, from 6 = 0.7 ± 0.2 for G = 2 to 6 = 0.03 ± 0.06 for 
G = 6 as for an integrable system (see section II of SM for 
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details). Thus, breaking the regularity of the fractal re¬ 
sults in the loss of the scale invariance of nearest-neighbor 
level spacings. 

As another control, we analyzed the equivalent regu¬ 
lar lattice of holes of equal size. This system retains the 
geometrical symmetries of the square but lacks the scale 
invariance since it has only two characteristic scales—the 
size of the holes and the size L of the system. We find 
that the spectrum of the lattice exhibits none of the re¬ 
markable properties of the fractal (scaling of the energy 
levels, presence of a regime with band-gap structure, self¬ 
similarity). This system must be desymmetrized in or¬ 
der to fit the spacings distribution to the form given by 
Eq. (1); the fitted distribution has a positive Brody pa¬ 
rameter that evolves to b ss 0 for the largest generations 
(see SM section II). Overall, these two negative controls 
suggest that the unique features of the spectrum, of the 
fractal are due to the combined effect of geometrical sym¬ 
metries and scale invariance. 

The electron gas. We used 1000 single-particle states in 
order to study the thermodynamic properties of the non- 
correlated electron gas [16] in a Sierpinski carpet. The 
Fermi level of the system for different total number of 
electrons N e with spin degeneracy is presented in Fig. 4. 
It is shown that the generation G can be chosen to ob¬ 
tain a Fermi level located close to a gap, and reproduce 
the properties of a semiconductor material with tailored 
band gap. The heat capacity C(T ) for N e = 100 elec¬ 
trons is shown in Fig. 4 (see section III of SM for details). 
The 1000 states suffice to describe each system accurately 
in a window of temperatures of almost three times the 
Fermi temperature Tp (Tp = nNp « 314 for G = 0). 
For G = 0 we recover the behavior of the 2D electron gas 
(2DEG), as expected. Indeed, C(T) depends linearly on 
T at low temperatures and converges to the classical limit 
( C c i = N e in two dimensions) at T ~ 2 Tp. For G = 1 
and G = 2 we obtain a similar behavior, with small de¬ 
viations from the behavior of the 2DEG. However, for 
G > 3 the presence of quasi-continuous bands separated 
by large energy gaps becomes important at low temper¬ 
atures, where the increasing temperature permits some 
electrons to suddenly cross the energy gaps and popu¬ 
late higher energy bands in a process that is reflected 
in the oscillations of the heat capacity as a function of 
temperature. For the highest generations, oscillations be¬ 
come substantial in a wider range of temperatures and, at 
higher temperatures, the heat capacity eventually over¬ 
shoots the classical limit. Remarkably, the dependence 
of the heat capacity per particle on the temperature ex¬ 
pressed in units of Fermi temperature converges to the 
same limit curve for the highest generations (Fig. 4). 

Two different mechanisms contribute to the overshoot¬ 
ing: On one hand, the observed oscillations due to the 
thermal population of high energy bands in regime (i) be¬ 
come more relevant as the generation increases. Figure 4 
shows that for N e = 100 at a temperature of T ~ 2 Tp 
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Figure 4. Top panel: Fermi energy Ep with respect the 
ground state energy of noninteracting electron gas for N e = 
100, 200, 500, and 1000 electrons for all generations studied 
here. Bottom panel: Dependence of the heat capacity per par¬ 
ticle C/N e on temperature (in units of Fermi Temperature) 
T/Tp for N e = 100 electrons. The alternation of bands and 
gaps generates oscillations in C as electrons thermally popu¬ 
late higher energy bands, and results in an overshooting with 
respect to the asymptotic thermodynamical value C = JV e 
for a 2DEG (shown as a horizontal black line), even at tem¬ 
peratures where the equivalent non-fractal 2DEG has reached 
the classical limit. Remarkably, it converges to a well-defined 
function of T/Tp as the generation increases. 

where the 2DEG (G = 0) already behaves classically, an 
analogous system of the same size but with the shape of 
a Sierpinski carpet still behaves quantum-mechanically. 
On the other hand, the increasing value of the ratio of 
the perimeter and surface of the fractal with increasing 
G affects the energy dependence of the density of states 
in regime (ii). Indeed, the first correction term to Weyl’s 
law is proportional to the perimeter and its dependence 
on \fe gives rise to the overshooting (see SM section IV), 
this effect becoming more pronounced as the perimeter 
term dominates over the surface term. However, for any 
finite value of G, the surface contribution eventually dom¬ 
inates over the perimeter correction at sufficiently high 
energies, and thus the asymptotic value agrees with the 
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2D ideal-gas value. This would not happen in the limit 
G —> oo (for a fixed energy or temperature); in this limit 
the surface of the fractal becomes zero while the perime¬ 
ter becomes infinite, yielding a non-classical result of a 
divergent heat capacity. 

In summary, the energy spectrum of quantum particles 
trapped in a Sierpinski structure shows self-similarity, 
with level attraction and repulsion leading to bands and 
gaps at any scale. The distribution of consecutive energy 
level spacings exhibits a transition from an exponential 
to inverse power law density as the fractal generation in¬ 
creases. At high generations, no characteristic size for the 
spacings can be defined and level spacings of any order 
of magnitude are equally probable. In fractal geometries, 
quantum effects originating from the discreteness of the 
energy levels can be observed at temperatures where the 
equivalent non-fractal system already behaves classically. 

In view of our results, fractal-shaped wires could ex¬ 
hibit promising features for a wide range of applications. 
Scale invariance of the spectrum could improve the pho¬ 
toelectric interaction of nanodevices in a wide range of 
light frequencies and allow the design of nanomaterials 
with tailored Fermi levels, band gaps, and conduction 
bands that can be thermally populated—mimicking the 
properties of semiconductors, as seen in the oscillations of 
the heat capacity—with potential applications for solar 
cells [30, 31], artificial photosynthesis [32], or metamate¬ 
rials [4, 5, 13, 33]. In addition, the persistence of quantum 
effects at high temperatures could increase the operating 
temperature range of these quantum nanodevices. 
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